Dynamic evolution characteristics and driving factors of tourism ecosystem health in China

Tourism ecosystem health is key to high-quality tourism development. China is now promoting sustainable development and high-quality transformation and upgrading of regional tourism; thus, the research on tourism ecosystem health is of practical significance. Based on the DPSIR model, an evaluation index system of tourism ecosystem health in China was constructed. Then the entropy weight method, spatial autocorrelation analysis, Markov chain analysis, and quantile regression were used to explore the dynamic evolution characteristics and driving factors of tourism ecosystem health in China from 2011 to 2020. The following conclusions were drawn: (1) The tourism ecosystem health in China showed an M-shaped fluctuation process as a whole, with significant spatial correlation and spatial difference. (2) There was a “path-dependent” and “self-locking” effect on the type transfer of tourism ecosystem health, and the type transfer was mainly between adjacent types in successive transfers, with the probability of downward transfer higher than upward transfer, and the geospatial background played a significant role in its dynamic evolution process. (3) In provinces with low tourism ecosystem health type, the negative effect of technological innovation capacity was more significant, and the influence coefficient of the positive effect of tourism environmental regulation and information technology level was larger, while in provinces with high tourism ecosystem health type, the negative effect of tourism industry agglomeration was more significant, and the influence coefficient of the positive effect of tourism industry structure and tourism land-use scale was larger.


. Introduction
The tourism business, like other industries, has a paradoxical connection with the ecological environment due to the dual industrial qualities of environmental dependence and resource consumption (1)(2)(3). With the rapid growth of tourism and the advent of "mass tourism", health issues affecting the industry's ecosystem have gradually surfaced (4,5). Phenomena like excessive resource consumption, ecosystem degradation, biodiversity loss, and rising environmental pollution limited the resilience of sustainable tourism development (6)(7)(8). It has become a paramount practical concern to find a solution to integrate the interaction between ecological protection and tourism development and sustain the health of the tourism ecosystem. Exploring sustainable tourism development models has emerged as the most crucial challenge in China's tourism development as one of the nations with the fastest global tourist growth rates (9,10). The report of the 20th CPC National Congress and the 14th Five-Year Plan for Tourism Development both place "promoting harmonious coexistence between man and nature" at the top of the overall agenda under the new model, which conforms to the normal distribution conditional mean, to identify the factors influencing tourism ecosystem health. Consequently, the possibility of different influencing factors in different regions of the tourism ecosystem health level was ignored. In addition, given the higher strategic value of sustainable tourism development and high-quality transformation and upgrading, little attention has been given to the spatial correlation and dynamic transfer of tourism ecosystem health at the national scale in the results of the available study, which are mostly focused on the dimensions of specific areas and economic zones. This paper attempts to deepen the previous research further based on these three deficiencies. Therefore, we selected 30 provinces (excluding Tibet, Hong Kong, Macao, and Taiwan) in China as the research subjects, then constructed a tourism ecosystem health evaluation index system based on the DPSIR model and used the entropy method to assess that health from 2011 to 2020. Meanwhile, the dynamic evolution characteristics of the tourism ecosystem health were explored through the entropy weight method, spatial autocorrelation analysis, and Markov chain analysis. Moreover, the panel quantile regression model was used to identify the driving forces behind the shifting trends in tourism ecosystem health under various quantile conditions. The research findings of this paper are intended to comprehensively understand the dynamic evolution characteristics and driving factors of tourism ecosystem health in China and provide a theoretical basis and decision-making reference for provinces with different levels of tourism ecosystem health so that it can promote the coordinated development of tourism and the eco-environment.
. Materials and methods

. . Entropy weight method
The subjective assignment approach introduces the impact of human factors. However, the entropy weighting method uses the original information of the indicators as the foundation for assigning weights, which can adequately reflect the significance of each indicator in the comprehensive index. The entropy weight method calculated the evaluation index weight of tourism ecosystem health. The following are the specific steps for implementation (35): The first step is to standardize the evaluation indexes.
Negative index : x ′ ij = max x ij − x ij max x ij − min x ij (2) In Formulas (1) and (2), x ij represents the original index j in the region i; x ′ ij represents the index after standardization processing; maxx ij and minx ij represent the maximum and the minimum values of index j, respectively.
The second step is to calculate the weight. In Formulas (3) and (4), w ij is the proportion of index j under the region i; e j is the entropy value of index j, and φ j is the index weight. The specific Frontiers in Public Health frontiersin.org . /fpubh. . Formula is as follows: The third step is a comprehensive evaluation index of tourism ecosystem health (TEH i ).
Integrating the results of the actual measurement of tourism ecosystem health in China and the synthesis of existing studies (28,33), the tourism ecosystem health of 30 provinces in China was divided into four types: non-healthy level, sub-healthy level, generally healthy level, and very healthy level (Table 1).

. . Global spatial autocorrelation analysis
Global spatial autocorrelation analysis, frequently measured by the global Moran's I index, is used to evaluate the overall spatial dependency between different geographical areas (36). In this paper, the global Moran's I index can reflect the spatial autocorrelation of tourism ecosystem health in China as a whole. The Formula for calculating the global Moran's I index is: In Formula (8), n represents the total number of provincial samples, which is 30 here; x i and x j represent the tourism ecosystem health measures of provinces i and j, respectively; x is the arithmetic means of the tourism ecosystem health of all provinces; W ij is the adjacency weight matrix, which indicates the adjacency relationship between two provinces, and W ij = 1 when provinces i and j are adjacent, otherwise it is 0. The value of global Moran's I index is [−1, 1]. The more the value of this index tends to 1, the stronger the correlation in the tourism ecosystem health space. The value of this index tends to 0, which indicates spatial decorrelation and random distribution in space.

. . Markov chain analysis method
The Markov chain, a transition matrix and the simplest stochastic model has been extensively used for state change studies at various spatial scales (37). In this paper, the Markov chain analysis method was used to calculate the initial and transfer probability of different states of tourism ecosystem health, and the trend over time was determined. Assuming that M t = [ M 1,t , M 2,t , M 3,t ,. . . M k,t ] is a vector of 1 × k state probability distributions in year t, the transfer between various types of tourism ecosystem health in different years can be represented by a Markovian transfer probability matrix of order z × z: In Formula (9), p ij represents the probability of a random process transitioning from state type i in year t to state type j in year t + 1, and it can be calculated as α ij /α i ; α ij represents the number of provinces of type i in year t transferred to type j in year t + 1 and α i represents the number of provinces of type i in all years.
To examine the relationship between the probability of transfer of tourism ecosystem health types and neighboring provinces, the spatial Markov chain analysis incorporated the spatial lag factors into the analytical framework based on the traditional Markov chain analysis (38). The traditional z-order transition probability matrix (z × z) is broken and decomposed into n conditional transition probability matrices (z × z × z) by the spatial Markov chain transition probability matrix. In the z-th condition matrix, P zij designates the probability that a province a that was in type i in year t transfers into type j in year t + 1, conditional on a spatial lag type of z. The spatial lag value for province a is a weighted average of the tourism ecosystem health of the province's spatial neighbors. The calculation Formula is: In Formula (10), W ab is the adjacency weight matrix, which indicates the adjacency relationship between two provinces, and W ab = 1 when provinces a and b are adjacent; otherwise, it is 0. Y b represents the tourism ecosystem health in province b; Lag a represents the spatial lag value of province a, indicating the state of the neighborhood of province a. By comparing the traditional Markov chain matrix and the spatial Markov chain matrix, it is possible to determine the importance of the surrounding area on the probability of a particular spatial province's type transition.

. . Panel quantile regression model
Koenker and Basett proposed the quantile regression model to overcome the shortcomings of the traditional regression method, which can only obtain the average influence of the explanatory variables on the dependent variable (39). The difference in significance levels and regression coefficients at various quantiles in the quantile regression model can reflect the heterogeneity of the influence of the explanatory variables on the explained variables. Tourism ecosystem health varies from province to province, and explanatory variables affect those with high tourism ecosystem health differently than those with low tourism ecosystem health. Therefore, the quantile regression model can comprehensively, systematically, and dynamically reveal the effect of driving factors on tourism ecosystem health. For the panel data, the quantile regression model is set as follows (40): In Formula (11) and (12), x it represents the independent variable of province i in the 1×k dimension of year t; β i and α i are the parameters to be estimated, and µ is a random error; Q yit represents the tourism ecosystem health at the quantile under the given explanatory variable conditions; τ represents the quantile points. The following Formula generally calculates the parameter β: In Formula (13), ρ τ represents the quantile loss function; W k represents the weight coefficient of the k quantile; β(τ q ) represents the influence coefficient of the k quantile.
Based on the current state of China's tourism industry and data availability, six indicators were selected to quantify the factors driving the change process in the spatial and temporal patterns of tourism ecosystem health under multiple factors. These indicators include tourism industry structure, tourism industry agglomeration, tourism environmental regulation, information technology level, technological innovation capacity, and tourism land-use scale ( Table 2).

. . Evaluation index system and data source
The PSR model served as the foundation for the DPSIR model, which was first proposed and utilized by the European Environment Agency (EEA) to provide a more thorough understanding of the interactions and feedback mechanisms between human activities and the biological environment (41). The DPSIR model for tourism ecosystem operations is as follows: for an extensive period, economic and social development and tourism demand act as driving forces (D) on the ecosystem, resulting in several pressures (P) on the environment and changing the health state (S) of the tourism ecosystem, which in turn has a variety of impacts (I) on individuals, nature, and society, forcing human society to respond to ecological changes (R). Combining the characteristics of tourism ecosystems and the concept of tourism ecosystem health, this paper constructs a set of indicators for assessing the tourism ecosystem health in China (13,42), as shown in Table 3. The data in this paper were mainly derived from the "China Statistical Yearbook", "China Tourism Statistical Yearbook", "China Culture and Tourism Yearbook", the yearbooks of each province, and the Statistical Bulletin of National Economic and Social Development of each province. Some missing data was supplemented using the linear interpolation method. Following data collection, panel data for 30 Chinese provinces from 2011 to 2020 were obtained.

. Results
. . Spatial and temporal characteristics of tourism ecosystem health . . . Time-series evolutionary characteristics According to Formula (1)-(7), the tourism ecosystem health of 30 provinces in China from 2011 to 2020 was evaluated based on the constructed evaluation index system of tourism ecosystem health, and the trend of change was plotted ( Figure 1). During the study period, the tourism ecosystem health in China showed an Mshaped fluctuation process of "increasing-decreasing-increasingdecreasing", with a three-stage evolution. In the first stage   The kernel density of the Gauss kernel function, depicted by Matlab R2011b software in Figure 2, was estimated to define the time-varying process of the absolute differences in the tourism ecosystem health across provinces. Then, utilizing the distribution location, distribution trend, and distribution extensibility, the kernel density was estimated using the Gauss kernel function to illustrate the dynamic evolution of the tourism ecosystem health in China. The nuclear density distribution curve shifted first to the right, then to the left, then to the right, and finally to the left from 2011 to 2020. This indicates that the tourism ecosystem health in China has undergone the process of "increasingdecreasing-increasing-decreasing", which is in line with the trend analyzed above. The overall kernel density exhibited a right-trailing phenomenon concerning the extension of the distribution over time, indicating that the gap between China's tourism ecosystem health and the average has widened and that the types of provinces with high tourism ecosystem health levels were increasing faster Frontiers in Public Health frontiersin.org . /fpubh. .

FIGURE
The trend of tourism ecosystem health in China from to . while the types of provinces with low tourism ecosystem health levels were decreasing. From a morphological point of view, the overall kernel density curve did not show a decreasing peak or increasing width. However, the sample period was characterized by multiple peaks, implying a particular gradient of differences in the tourism ecosystem health in China.

. . . Spatial di erentiation characteristics
According to Formula (8), combined with the measured value of tourism ecosystem health, using the Rook spatial weight matrix (Hainan and Guangdong were set as neighbors to avoid the "island phenomenon"), the spatial autocorrelation of tourism ecosystem health in China is tested and analyzed, and the global Moran's I index test results were calculated using Stata 16.0 software (Table 4). It can be seen from Table 4 that the global Moran's I index of tourism ecosystem health in China from 2011 to 2020 was all positive and passed the statistical test, reflecting a solid spatial correlation. China's spatial distribution of tourism ecosystem health was not isolated and random. However, it showed apparent spatial dependence and agglomeration, such as convergence between provinces with high tourism ecosystem health level types, forming a "high-high" agglomeration pattern, and the convergence between provinces with low tourism ecosystem health level types, forming a "low-low" agglomeration pattern. From 2012 to 2020, Moran's I index of tourism ecosystem health continued to decline in fluctuation, showing a trend of weakening spatial dependence, which, to a certain extent, indicated that China's tourism ecosystem health was gradually showing a coordinated and linked trend of regional integration.
With the help of ArcGIS 10.2 software, the spatial distribution map of tourism ecosystem health in 2011, 2014, 2017, and 2020 was drawn, and the tourism ecosystem health was divided into four types, which were non-healthy level, sub-healthy level, generally healthy level and very healthy level from low to high (Figure 3). During the inspection period, the provinces with higher tourism ecosystem health levels were mainly concentrated in Beijing, Shanghai, Guangdong, Jiangsu, and Zhejiang, accounting for a small portion of the total. This showed significant regional differences in the tourism ecosystem health in China, which were out of balance. Specifically, 28 provinces with non-healthy tourism ecosystem health levels in 2011, except for Beijing and Guangdong. Compared with 2011, Jiangsu and Zhejiang joined the provinces with generally healthy tourism ecosystem health levels, while Liaoning, Shandong, Fujian, and other provinces withdrew from the non-healthy level. In 2017, the provinces of the very healthy level type remained unchanged. Tianjin and Fujian went up from the sub-healthy level to the generally healthy level. In contrast, the provinces of the sub-healthy level type further expanded, most of which were concentrated in the central and western regions. In 2020, the spatial pattern of tourism ecosystem health in China was similar to that in 2017, with only minor changes found in the type of some provinces. From the spatial distribution pattern, it can be seen that the tourism ecosystem health showed the distribution law of decreasing gradually from the eastern coastal regions to the central, western, and northeastern regions. The provinces with higher tourism ecosystem health levels, driven by both their location and policy advantages, have achieved high-quality growth in the tourism economy under the intensive mode while simultaneously optimizing the structure of the tourism industry and environmental management through advanced technological advantages. The provinces with lower tourism ecosystem health levels were mainly severely constrained by their tourism resource endowment and location conditions, which, together with their strong dependence on tourism resources, has led to bottlenecks in the development of the tourism economy.

. . . Dynamic evolutionary characteristics
This section used the Markov chain method to investigate the state transition of provinces in the distribution of tourism ecosystem health. It attempted to explore the long-term dynamic change trend of the differences in tourism ecosystem health between the various provinces in China. According to the type classification, the tourism ecosystem health of provinces can be discretized into four state-level types. Namely, non-healthy level (0, 0.30], sub-healthy level (0. 30  respectively. The definition of an upward transfer is from a low level to a high level, and the definition of a downward transfer is from a high level to a low level. To further understand the dynamic evolution trend of tourism ecosystem health in China, the Markov transfer matrix of tourism ecosystem health in China was obtained using Matlab R2011b software ( Table 5). The value on the diagonal line represents the possibility that the province will always remain at a particular statelevel type. In contrast, the value on the off-diagonal line represents the probability that the various provinces will transition from one state-level type to another. Therefore, without considering the spatial effect, the dynamic evolution of tourism ecosystem health is characterized as follows: (1) The probability value on the diagonal line in the Markov transition matrix was significantly larger than those on the non-diagonal, with a minimum value of 72.41% and a maximum value of 100.00%. Under the above conditions, tourism ecosystem health attained at least 72.41% probability in the future development process and remained at the same state-level type. The probability of maintaining stability at each state-level type was greater than the probability of the upward transfer and the downward transfer, which reflected that the club convergence of tourism ecosystem health types was highly .
/fpubh. . stable, and its type transfer had a "path-dependent" and "selflocking" effect. (2) In terms of the probability value on both sides of the diagonal line, the probability of the non-healthy level type transitioning upward was 27.58%. It can be seen that the provinces with non-healthy level types were more likely to transfer upward, indicating that the tourism ecosystem health in the nonhealthy level provinces was unstable and prone to state leapfrogging compared to other states. The probability of the sub-healthy level type transitioning downward was 11.88%, and the probability of the type transitioning upward was 7.92%. The probability of the downward shift of tourism ecosystem health at the sub-healthy level was more significant than the probability of the upward shift, which reflected that tourism ecosystem health still had a negative trend. The probability of the generally healthy level type transitioning downward was 27.59%, with the risk of a precipitous fall. (3) The probability values on both sides of the non-diagonal line were significantly smaller than those on both sides of the diagonal line, with a maximum value of only 0.86%, indicating that the improvement of tourism ecosystem health was a gradual process and that it was difficult to achieve leapfrogging in the short term. The issue of spatial correlation is disregarded in the traditional Markov chain approach because it presumes that regions are independent of each other. This paper incorporated spatial lag into the traditional Markov chain analysis to investigate the neighborhood's influence and the probability that it will change among convergent groups. The results of the spatial Markov chain transition probability matrix of tourism ecosystem health in China from 2011 to 2020 are presented in Table 6.
The issue of spatial correlation is disregarded in the traditional Markov chain approach because it presumes that regions are independent of each other. This paper incorporated spatial lag into the traditional Markov chain analysis to investigate the neighborhood's influence and the probability that it will change among convergent groups. The results of the spatial Markov chain transition probability matrix of tourism ecosystem health in China from 2011 to 2020 are presented in Table 6.
Through comparison with Tables 5, 6, the following spatial dynamic evolution characteristics of tourism ecosystem health could be obtained after considering the geospatial background: (1) Geospatial background played a significant role in the dynamic evolution process of tourism ecosystem health. The transfer probability of tourism ecosystem health varied against the neighborhood background, implying different levels of health. It also differed from that calculated using the traditional Markov probability transfer matrix. In particular, the probability of maintaining stability at the non-healthy level type in the traditional Frontiers in Public Health frontiersin.org . /fpubh. .
Markov probability transfer matrix was 72.14%. In comparison, the probability of maintaining its original state-level type on different geospatial backgrounds was 72.73, 72.50, 87.50, and 0.00%, respectively, indicating a significant difference between the probability transfer of tourism ecosystem health with and without considering the geospatial background.
(2) China's tourism ecosystem health exhibited a "spatial spillover" effect, with the tourism ecosystem health levels affecting each other between neighboring provinces. In general, the probability of the statelevel type of province transferring downward would increase if it was adjacent to a province with a low tourism ecosystem health level. In contrast, the probability of the state-level type of province transferring upward would increase if it was adjacent to a province with a high tourism ecosystem health level. Specifically, as the state-level type of neighborhood tourism ecosystem health rose, the probability of the provinces of the non-healthy level type transitioning upward was 27.28, 27.50, 12.50, and 100.00%, showing a fluctuating upward trend; the probability of provinces of the generally healthy level type transitioning downward was 100.00, 50.00, 20.83, and 0.00%, showing a precipitous downward trend. This suggests that provinces with higher tourism ecosystem health levels have a radiation effect on their neighbors, whereas provinces with lower tourism ecosystem health levels may have an inhibition effect. Therefore, provinces with a high tourism ecosystem health level need to exploit their spillover effect actively and, through their radiation capacity, promote the health type of the surrounding provinces by leading from point to point.

. . Analysis of driving factors
In order to ensure the smoothness of the panel data time series and reduce heteroscedasticity in this paper, natural logarithms were taken as all non-percentage variables. At the same time, this paper employed both the LLC test, the Breitung test (which assumes a common unit root process), the IPS test, the Fisher-ADF test, and the Fisher-PP test (which assumes an individual unit root process) to probe the unit root properties of the study variables so that it could avoid problems such as multicollinearity and spurious regressions. According to the estimation results, every variable strongly rejected the null hypothesis that the panel had a unit root at a 1% significance level, indicating that the panel data were smooth and could be suitable for further regression analysis. This paper selected five representative quantiles for analysis, including 10, 25, 50, 75, and 90%. The Markov chain Monte Carlo (MCMC) method was adopted to estimate the driving factors of tourism ecosystem health to avoid possible endogeneity. The estimation results were presented in Columns (1) to (6) of Table 7. The changes in the quantile regression coefficients at different quantiles were displayed graphically in Figure 3. The estimated results were from stata16.0 software.
As shown in Table 7 and Figure 4, as a whole, tourism industry agglomeration and technological innovation capacity imposed negative influences on tourism ecosystem health, and tourism industry structure, tourism environmental regulation, information technology level, and tourism land-use scale were all significant in influencing tourism ecosystem health. However, the influence coefficients of each variable at different quantiles were significantly different.
(1) The results of the quantile regression of tourism industry structure (lnTIS) showed an M-shaped pattern, which was always positive and significant, indicating that the impact of tourism industry structure on tourism ecosystem health varies at different quantiles. The government should therefore transform the tourism economy and promote the integration of tourism industries to enhance tourism ecosystem health. It was important to note that the coefficient of influence of the tourism industry structure on the high quantile was relatively large. Therefore, the government should pay more attention to and strengthen the upgrading of the tourism industry structure in provinces with high tourism ecosystem health, thereby reducing duplicate construction and resourceconsuming projects and increasing low-consumption, highquality service projects, thereby promoting sustainable tourism development in these provinces. (2) At the 10, 50, 75, and 90% quantiles, the absolute value of the coefficient of tourism industry agglomeration (lnTIA) had an increasing trend, indicating that tourism industry agglomeration was the main obstacle factor for provinces with high tourism ecosystem health in China, that is, the higher the level of tourism ecosystem health, the greater the pressure of tourism industry agglomeration on the ecosystem. Consequently, our results suggested that strengthening the centralization of resource use and the use of large-scale pollution control infrastructure in provinces with high tourism ecosystem health was one of the most effective ways to increase the resilience of sustainable tourism development in such provinces. (3) Tourism environmental regulation (TER) significantly affected tourism ecosystem health at the lower quantiles, but this had no significant marginal impact as the quantiles increased. This may be because of the spatial mismatch between supply and demand in China's tourism resources and natural ecological background. The provinces with higher tourism ecosystem health types are generally more market-oriented regions, where the market itself can achieve an efficient allocation of resources through competitive mechanisms, price mechanisms, and supply and demand mechanisms. Too much macro-regulation is not conducive to forming environmentaleconomic systems such as green capital markets, green credit, ecological compensation, and other environmentaleconomic systems. (4) The influence of technological innovation capacity (lnTIC) on tourism ecosystem health was negative, with a non-significant negative effect on technological innovation capacity at the higher quartiles (75% and 90%). All other quartiles showed a significant adverse effect at the 1% level. The absolute value of the coefficient of technological innovation capacity increased as the quantile decreased, with the negative effect reaching a maximum at the 10% quantile. According to our findings, the main obstacle factor for provinces with low tourism ecosystem health in China was a lack of technological innovation capacity.   t statistics in parentheses. * P < 0.05, * * P < 0.01, and * * * P < 0.001.

FIGURE
Changes in quantile regression coe cients for tourism ecosystem health.
coefficients for the first 50% of the quantile being significantly larger than those for the second 50% of the quantile. From the coefficients, the degree of impact of information infrastructure and information technology consumption was more substantial for provinces with low tourism ecosystem health. In contrast, for provinces with high tourism ecosystem health, the effect of increasing informatization was diminished. does not exceed the reasonable carrying capacity of the tourism environment.

. Conclusions and discussion
Based on the DPSIR model, we systematically constructed an evaluation index system to measure China's tourism ecosystem health. Then the dynamic evolution characteristics and driving factors of tourism ecosystem health in China from 2011 to 2020 were analyzed with the entropy weight method, spatial autocorrelation analysis, Markov chain analysis, and quantile regression. The main conclusions are as follows: (1) During the ten years from 2011 to 2020, the tourism ecosystem health in China showed an M-shaped fluctuation process as a whole, with three distinct stages of rapid increase (2011-2012), slow decrease (2012-2015), and fluctuating development (2015-2020). The kernel density exhibited a right-trailing phenomenon and multiple peaks, indicating that the gap between China's tourism ecosystem health and the average has widened and that there was a specific gradient of differences. Regarding spatial differentiation, the tourism ecosystem health in China was significantly spatially correlated, characterized by clustered and contiguous development, with significant regional differences in distribution. (2) There was a "path-dependent" and "selflocking" effect on the type transfer of the tourism ecosystem health in China. Type transfer generally occurs between adjacent types in successive transfers, while the probability of cross-type transfer is small, and the probability of downward transfer is higher than upward transfer. In addition, geospatial patterns played a significant role in the dynamic evolution process of tourism ecosystem health. Specifically, the probability of the state-level type of the province transferring downward would increase if it was adjacent to a province with a low tourism ecosystem health level, while the probability of the state-level type of the province transferring upward would increase if it was adjacent to a province with a high tourism ecosystem health level. (3) The tourism ecosystem health in China was driven by a combination of factors, including tourism industry structure, tourism industry agglomeration, tourism environmental regulation, information technology level, technological innovation capacity, and tourism land-use scale. Moreover, in provinces with low tourism ecosystem health type, the dampening effect of technological innovation capacity was more significant, and the influence coefficient of the positive marginal effect of tourism environmental regulation and information technology level was larger, while in provinces with high tourism ecosystem health type, the dampening effect of tourism industry agglomeration was more significant, and the influence coefficient of the positive marginal effect of tourism industry structure and tourism land-use scale was larger.
Under the new development concept of "innovation, coordination, green, openness and sharing", a systematic and in-depth study on the dynamic evolution of tourism ecosystem health in China and its driving factors are of great significance in promoting the sustainable development and high-quality transformation and upgrading of regional tourism. The main contributions of this paper are listed as follows: (1) The dynamic transfer process and pattern of tourism ecosystem health in each province of China from 2000 to 2015 were analyzed using the spatial Markov chain analysis, which can visually reveal the heterogeneity of the "spatial spillover" effect of tourism ecosystem health and the influence of geospatial background. (2) In contrast to the idealistic treatment model of mean regression, the panel quantile regression model emphasized the heterogeneity of driving factors in the context of various tourism ecosystem health types. Its empirical findings could more accurately reflect the actual situation. Thus, the results can provide a research methodological reference for a more comprehensive, systematic, and dynamic exploration of the driving mechanisms of tourism ecosystem health in similar areas, especially in developing countries, in the future. (3) The research scale was reduced to the provincial level, and the heterogeneity and regularity of tourism ecosystem health at the regional scale could be better explained, providing empirical support for local governments to formulate appropriate tourism ecosystem health strategies for different provinces.
The findings of this paper have important policy implications: (1) The government and tourism authorities should give due consideration to regional synergy and integrated management in tourism cooperation and development and tourism environmental protection policies as a means of reducing the constraints and impacts of spatial effects on tourism ecosystem health, reducing spatial differences in tourism ecosystem health between provinces, and achieving the goal of regional coordination and sustainable development of the tourism industry.
(2) More attention should be paid to the dynamic evolution of tourism ecosystem health in different provinces, especially in provinces with a "generally healthy" neighborhood type, to avoid the risk of downward transfer resulting from the crude growth of their tourism economies. (3) The provinces with high tourism ecosystem health must make efforts to accelerate the process of allocating tourism industry elements, continuously promote structural reform on the supply side of tourism, guide tourism development and transformation and upgrading with the concept of ecological priority and green development, and actively cultivate new low-carbon and green tourism industries. Meanwhile, actions must be taken to initiate the allocation of land resources in ways that are compatible with the direction that the tourism industry is going as well as through intensive utilization for sustainable development. For provinces with low tourism ecosystem health, the research and development of tourism pollution treatment and prevention technologies should be focused on further improving tourism ecosystem health. Additionally, regional tourism ecosystem health policies should be continuously improved. Tourism ecosystem health should be promoted through new media platforms to encourage green travel and low-carbon consumption among tourists.

Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.